Finite Size Scaling Behavior of Dissipative Tunnel Junctions 
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We present the results of a series of Quantum Monte Carlo calculations of the temperature depen- 
dent conductance of dissipative non-superconducting tunnel junctions. Finite size scaling methods 
are used to demonstrate the absence of coherent transport in Ohmic tunnel junctions. However, 
a quantum phase transition between coherent and incoherent behavior is found in sub-Ohmic tun- 
nel junctions. The critical conductance is found to be temperature independent and the critical 
exponents, and r), are found to be in agreement with renormalization group results by Kosterlitz. 
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In spite of the significant effort to understand the effect 
of ohmic dissipation on normal (non-superconducting) 
tunnel junctions, there remains a conflict between an- 
alytical studies and quantum Monte Carlo simu- 
lations |^,|| of these systems. In particular, analytical 
treatments such as renormalization group [H, large N 
approximation ^ , and variational methods ^ generally 
conclude that no quantum phase transition occurs in tun- 
nel junctions subjected to a purely Ohmic environment. 
In contrast, it has been argued Q that quantum Monte 
Carlo studies of Ohmic tunnel junctions exhibit a quan- 
tum phase transition between an ordered phase which ex- 
hibits coherent tunneling and a disordered phase which 
does not. 

If such a phase transition existed, it would be a remark- 
able example of a quantum phase transition driven by 
dissipation. Moreover, the existence of such a transition 
would be of considerable importance to the electronic 
and transport properties of a variety of systems including 
double quantum dots and metallic granular composites. 
For these reasons, we believe that it is essential to resolve 
this discrepancy. Towards this end, we will reexamine the 
simulation data using using finite size scaling methods. 
This approach will allow us to identify quantum phase 
transitions which may occur in Ohmic and sub-Ohmic 
tunnel junctions. 

Our discussion is based on the generalized Ben Jacob 
Mottola Schon (BMS) model which describes the ef- 
fects of dissipation and Coulomb blockade phenonenon on 
non-superconducting tunnel junctions. The BMS model 
is defined by the effective action 



dr dr'air - t')[1 - cos(<?!)(t) - 0(t'))] (1) 



where the phase 4){t) is defined in terms of the time- 
dependent voltage across the tunnel junction, V{t), such 
that V{t) = f0(T). The generalized BMS model is a 
long range XY model with spin-spin interactions of the 

form a(T) = aoTg' [ ^i^^J^Tr) ) ■ defined in 

ey2C by TQ^n/EQ. 



^ sin(7rkBTT) ^ 

terms of the charging energy Eq - 
Charging effects are introduced to the model by working 



on a lattice such that 0(r) is defined for r = Jtq where 
j = 1...N and where N = Ti(3/tq = Eq /ksT is the total 
number of timeslices. e is a parameter, discussed below, 
which characterises dissipative bath coupled to the tunnel 
junction. When e = 0, the coupling parameter, a^, is 
related to Rt, the T ^ oo limit of the resistance, by 
the expression ag = ^-/(ZTre^i?^). More generally ag 
is proportional to \t\'^NLNji i.e. the squared tunneling 
matrix element multiplied by the density of states of the 
left and right electrodes. 

When one derives this action from a more microscopic 
model, the long range interaction, a(T), is generated 
by intergrating out the modes associated with a dissi- 
pative bath. These modes might, for instance, repre- 
sent, particle-hole excitations occuring in the tunnel junc- 
tion electrodes. In the absence of manybody effects, the 
particle-hole excitations give rise to Ohmic dissipation 
(e = 0). However, in some situations, Fermi edge singu- 
larity effects such as the orthogonality catastrophe and 
the exciton effect P,[lO|| may become important. [pT|J^. In 
such cases, the dissipative bath is not Ohmic i.e. e 7^ 
and the tunnel junction may exhibit a variety of interest- 
ing effects including non-linear I{V) characteristics. 

As discussed by Drewes et al |ll|], the renormaliza- 
tion group theory predicts a quantum phase transition 
for sub-Ohmic tunnel junctions (i.e. e > 0). This 
transition occurs at ao = ade) where ac(e) ^ cxd as 
e ^ 0. This result appears to be in conflict with the 
Monte Carlo simulations of Ohmic (e = 0) tunnel junc- 
tions by Scalia et al The latter study interpreted 
the observed ao dependence of the correlation function 
(7(r) =< expi(/)(T) exp — i(/)(0) > as evidence for a quan- 
tum phase transition at « 0.7. 

In this paper, we will reconsider the simulation results 
using finite size scaling methods. Using these methods 
we have obtained the following results: (1.) No quantum 
phase transition occurs in Ohmic tunnel junctions with 
ao < 10.0; (2.) The sub-Ohmic tunnel junctions exhibit 
a quantum phase transition characterized by a tempera- 
ture independent critical conductance of order e^//i; (3.) 
Values for the 77 and v critical exponents are found to 
be consistent with the analytical results of Kosterlitz 0] 
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and Nickel et al 



A Quantum Phase Transition: We begin our discussion, 
by considering the tunnel junction conductance which 
can be written as [|l]],f2l 



G 



27rao 



h(3 



7.(t)(cos(</.(t)-0(O))) dr (2) 



where 7,(r) = [nikBT/EQ)/ sm{nkBTTh)]-'. We evalu- 
ate G using gir) = (cos((/)(t) — 0(t'))) obtained from a 
series of simulations. The simulations involved the stan- 
dard Metropolis Rosenbluth Teller Teller algorithm of the 
long range XY model on a 1-D lattice of = TiP/tq spins 
and periodic boundary conditions. Our simulations in- 
cluded lO'^ cycle runs for the 10, 20, and 32 timeslice 
systems and 2 x 10^ cycle runs for the 64 and 128 times- 
lice systems. 

The results for e — 0.2 are presented in fig. |l| and the 
results for e = are presented in fig. ||. Interestingly, 
in fig. the N dependence of the conductance curves 
reverses at the point of intersection at Uc ~ 0.9. For 
Q!o > Oic one observes metallic behavior i.e. a regime 
where G increases as ksT — Eq/N — > 0. Conversely, 
for ao < Oic, G decreases as T ^ 0. Hence, we may 
identify the crossing point as a quantum phase transition 
between two phases. Hereafter, these two phases will 
be referred to as sub-Ohmic ( ao > etc) or insulating 
(ao < oic). The sub-Ohmic phase exhibits a conductance 
which increases with increasing N oc l/ksT whereas the 
insulating exhibits a conductance which decreases with 
increasing N. 

Remarkably the conductance curves in Fig. don't 
simply exhibit a set of crossings. Instead all six curves 
cross at a single point. This observation suggests that the 
critical conductance is temperature independent. The 
existence of a temperature independent critical conduc- 
tance can be understood using the finite size scaling 
ansatz 



< e 



where d = 1 is the dimensionality of the model, tq = 
h/Eq is the size of the time slices, and where is the 
correlation time. Eq. ^ together with eq. || determine 
the scaling behavior of the conductance: 



knT 



(4) 



where A = ?i/^t and where F+{x) and F-{x) are the 
ao < Oic and ao > olc branches of some universal scaling 
function. Now Fisher, Ma and Nickel have shown that 
the long range spin-spin interactions cause rj = 1 -\- e. 
Because of this, eq. ^ implies a temperature independent 
critical conductance, Gc{e), of order e'^/h. This result 
has been confirmed using a large N approximation of the 
model H]. According to that treatment 
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FIG. 1. Conductance for e = 0.2. Notice that the conduc- 
tance curves cross at Gc ~ 2le^/h. The crossing separates 
a phase where the conductance exhibits sub-Ohmic from a 
phase where it exhibits insulating behavior. 
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FIG. 2. Conductance for e = 0.0. 



G.(e)=2^(l-6)ctn(y)^ 



(5) 



This result shows that Gc —> as e 1 which indicates 
that the disordered phase is absent at e — 1. Similarly 
eq. H shows that Gc ^ cxo as e ^ 0. This suggests that 
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the ordered phase is absent at e = 0. These interpre- 
tations are supported by the observed e dependence of 
ac(e) given in fig. |. 

As was mentioned in the introduction, the absence of 
an ordered phase in the e — > limit is not only a predic- 
tion of the large N approximation but also is predicted by 
Kosterlitz' renormalization group treatment of the long- 
range XY model. ^ Additional evidence may be found 
in figs. 1^ and ^. For 10 > ao > 0.1, G decreases with 
decreasing temperature, ksT = Eq/N. This strongly 
suggests that no metallic or coherent phase is occurs in 
this range of ao. In addition, these figures do not exhibit 
any crossings which might indicate a phase transition. 
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FIG. 3. Critical value of ao vs e as determined by the cross- 
ing of conductance curves. 



Finite Size Scaling Behavior: We have seen that the ex- 
istence of the temperature independent conductance pro- 
vides indirect evidence of the finite size scaling behavior 
proposed in eq. ^. To obtain more direct evidence for 
finite size scaling, one must be able to to collapse the 
conductance data using 



G = 



h 



F± (^N^/" 



ao 



(6) 



This is done in fig. |[ In that figure, we used the value of 
ac from fig. |l| and adjusted v to minimize the scatter. In 
this manner we obtained data collapse with ly = 4.2 ± 0.6 
for e = 0.2 This may be compared with the one loop 
renormalization group result that ly = l/e[l -I- 0{e)]. 

Next consider the e = limit. In this limit, the scaling 
behavior indicated in cqn. a reduces to 
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FIG. 4. Same data as in previous figure, except G is plot- 
ted vs l/ao- Observe that the conductance is monotonically 
increasing function of temperature , ksT oc Eq/N. 



where ^ = lim^^o otc/v is assumed to be finite. The re- 
sulting data collapse occurs is displayed in fig. ^. A 
comparison between figs. ^ and ^ indicates that only the 
monodecreasing branch of the scaling function is present 
in the e = data. This means that for all a^, the conduc- 
tance of Ohmic tunnel junctions decreases with increas- 
ing N oc l/ZcsT. Again we conclude that, as predicted by 
Kosterlitz's renormalization group treatment, a coherent 
(ordered) phase is does not occur in Ohmic tunnel junc- 
tions. 

Summary: We have performed a finite size scaling anal- 
ysis of the conductance of disspative Ohmic and sub- 
Ohmic tunnel junctions. Consistent with the predictions 
of renormalization group and large N theories, the sub- 
Ohmic tunnel junction exhibits a quantum phase transi- 
tion between sub-Ohmic and insulating phases. The re- 
sults for Ohmic tunnel junctions are also consistent with 
the RNG predictions: In particular, little evidence for 
a metallic phase is found in Ohmic tunnel junctions for 
Qfo < 10. However, the conductance does exhibit finite 
size scaling consistent with a quantum phase transition 
which, according to renormalization group theory and 
large N approximation, occurs at a^T^ = in the Ohmic 
(e = 0) limit. 
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FIG. 5. The conductance ai e — 0.2. Two branches occur. 
The upper branch is metaUic (increasing with N) and the 
lower branch is insulating (decreasing with N). Data collapse 
vahdates scaling form in eq. H. 
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FIG. 6. Conductance for e — 0.0. Note: The conduc- 
tance is monotonicaUy increasing function of temperature, 
ksT — Eq/N. This implies insulating behavior at low tem- 
peratures. 
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